Relation between directed polymers in random media and random bond dimer models 
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We reassess the relation between classical lattice dimer models and the continuum elastic descrip- 
tion of a lattice of fluctuating polymers. In the absence of randomness we determine the density and 
line tension of the polymers in terms of the bond weights of hard-core dimers on the square and the 
hexagonal lattice. For the latter, we demonstrate the equivalence of the canonical ensemble for the 
dimer model and the grand-canonical description for polymers by performing explicitly the contin- 
uum limit. Using this equivalence for the random bond dimer model on a square lattice, we resolve 
a previously observed discrepancy between numerical results for the random dimer model and a 
replica approach for polymers in random media. Further potential applications of the equivalence 
are briefly discussed. 
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I. INTRODUCTION 

Dimer coverings of different lattice types have been 
employed recently as a starting point to study more com- 
plex physical systems like quantum dimer modelsi^, geo- 
metrically frustrated Ising magnets with simple quantum 
dynamics induced by a transverse magnetic field^ and 
elastic strings pinned by quenched disordes^. The com- 
mon concept of these approaches is to add to a classical 
dimer model with a hard-core interaction a perturbation 
in form of simple quantum dynamics, quenched disorder 
(random bonds) or additional (classical) dimer interac- 
tions. The classical hard-core dimer model can be solved 
on arbitrary planar graphs^'* . For bipartite lattices there 
exists a representation of the dimer model in terms of a 
height profile of a two-dimensional surfac o^d" . Steps sep- 
arating terraces of equal height form a lattice of directed 
and non-crossing polymers^. This polymer lattice has 
been used to study the effect of quantum fluctuations for 
a Ising antiferromagnet on a triangular lattice in a trans- 
verse field^*^*ii. Moreover, the dimer model with random 
bond energies can be simulated in polynomial time and 
provides an independent test of the replica theory for 
the pinning of elastic polymers by quenched disorder—. 
Hence it is important to understand the relation between 
lattice dimer models and the elastic continuum descrip- 
tion of the corresponding polymers and to relate the pa- 
rameters of the two models. 

Here we show that the canonical dimer ensemble in the 
classical case maps to the grand canonical ensemble of the 
polymer system. Regarding the polymers as imaginary- 
time world lines of free fermions in one dimension, we 
give a simple derivation of the continuum free energy of 
lattice dimer models which agrees with the exact result if 
the continuum limit is taken properly. Applying the re- 
lation between dimers and polymers to the dimer model 
with random bond energies, we resolve a previously ob- 
served discrepancy between numerical simulations of the 
dimer system and a replica theory for the polymers. We 



analyze which quantities of the pinned polymers can be 
probed by simulations of the random dimer model. More 
specifically, we study the dimer model both on the hexag- 
onal and the square lattice, and discuss the meaning of 
the different lattice symmetries for the polymer repre- 
sentation. In the presence of bond disorder, we focus 
on the square lattice since its polymer density is con- 
served independently of the disorder configuration, and 
hence shows no sample-to-sample variations. Our re- 
sults should provide a starting point for other situations 
where no exact solution of the dimer model is possible as, 
e.g., in a recently studied case of nearest neighbor dimer 
interactions^^. The analogy between directed polymers 
in two dimensions and Luttinger liquids could be applied 
to understand more general interacting dimer models. 

II. MODELS: DIMERS AND POLYMERS 

The partition function of the dimer model on the 
hexagonal lattice is given by 

zo = ^zr4% (1) 

{D} 

where the sum runs over all complete hard-core dimer 
coverings of the lattice, and ni and 712 are the numbers 
of dimers occupying the two types of non-vertical bonds 
of the hexagonal lattice, see Fig^ which carry weights zi 
and Z2, respectively. The weights on the vertical bonds 
are assumed to be unity. Using the same notation, we 
define for the square lattice the partition function as 

^□ = E^^' (2) 

{D} 

where ni is now the number of dimers covering horizontal 
bonds which have a weight of zi while all vertical bonds 
carry unity as weight. For these clean dimer models the 
partition function and correlation functions are known 
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from exact results^i^. For random weights, only numer- 
ical results are available, see, e.g., Refs. 00. However, 
there is a useful connection between the dimer models 
and non-crossing directed polymers in (1-1-1) dimensions. 
This relation is independent of the actual bond energies, 
and hence applies also to random dimcr models. This is 
particularly interesting since the random bond energies 
translate to a pinning potential for the polymers, a prob- 
lem whose continuum version can be studied in (1-1-1) 
dimensions by a replica Bethe ansatai^ii^. Numerical al- 
gorithms for the random dimer model provide hence a 
unique opportunity to probe the replica symmetric the- 
ory which is commonly used to describe pinning of elastic 
media. 

The relation between dimers and polymers is estab- 
lished by superposing every dimer configuration by a 
fixed reference dimer configuration. For the hexagonal 
lattice, the reference state consists of a covering of all 
vertical bonds, whereas for the square lattice a staggered 
covering of the vertical bonds is chosen, see Fig. ^ In 
the superposition state, a bond is covered by a dimer 
(or polymer segment) if either it is covered only in the 
original state or only in the reference state. Because of 
the hard core constraints for the dimers, the polymers 
are non-crossing and they are oriented along the vertical 
direction due to the choice of the reference state. 

The resulting lattice of polymers can be described in 
the continuum limit by an elastic theory which is of the 
form 
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with compression modulus cn, tilt modulus C44 and local 
polymer density p{r) — J2j ^{^ ~ Xj{y)), where Xj{y) is 
the path of the j^^ polymer. The random bond energies 
are accounted for by a random pinning potential ^(r) 
which is uncorrelated, i.e.. 



V{r)V{r') = AS{r - r' 



(4) 



so that A measures the strength of disorder. In order to 
compare results for the dimer and the polymer model, we 
establish a relation between the dimer weights and the 
elastic constants and mean density of the polymers. Let 
us consider first the clean limit with A = 0. It is obvious 
from the mapping between dimers and polymers that the 
polymer density can vary with the dimer covering. For 
example, if the dimer state matches exactly the reference 
state, the polymer density is zero. However, one can 
define a mean polymer density by averaging over all dimer 
coverings. 

For the hexagonal lattice the mean density is deter- 
mined by the mean number of occupied non-vertical 
bonds in the original dimer configuration so that po = 
(ni -|- 712) / (V^bQN) where N is the total number of 
dimers, Bq the lattice constant and (. . .) denotes here 
an ensemble average over all dimer coverings. This 
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if 2i -|- 22 > 1 and pQ = if 21 -I- 22 < 1. For the square 
lattice, the mean number of polymers is determined by 
the probability that a vertical bond is occupied by a seg- 
ment of the polymer. This is case if the bond is cov- 
ered by a dimer in the original dimer configuration and 
not covered in the reference state, or in the other way 
around. The probability that a vertical bond is covered 
by a dimer in the original covering is pd — 1/2 — ^(21) 
with (f>{zi) — arctan(2i)/7r. For the reference state it is 
simply = 1/2. Hence after the superposition of the two 
dimer states, the probability that a vertical bond is cov- 
ered by a polymer is given by j3rf(l—pr)+Pr(l—Pd) = 1/2 
independent of 21. This fixes the mean density ali^ 



Pn 
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where bn is the lattice constant. Notice that the den- 
sity on the square lattice does not change with the bond 
weight zi. 

The elastic constants are length scale dependent due 
to renormalization effects from the non-crossing con- 
straint. Since their is no additional interaction between 
the dimers than the hard core repulsion, the compres- 
sion modulus cii is zero on microscopic scales. A finite 
macroscopic cn is generated by a reduction of entropy 
due to polymer collisions, see below. The tilt modulus 
C44 — gp on microscopic scales (or at very low density) 
is given by the line tension 5 of a single polymer and the 
mean polymer density. The reduced line tension g/T of 
an individual polymer at temperature T can be obtained 
from a simple random walk on the lattice^ which per- 
forms transverse steps according to the weights of the 
dimer model. For the hexagonal lattice it readsii 



go ^ 2 + r, + l/Tj 
T 2bo 



with ri = 21/22, and for the square lattice one has 

5n _ 221 + 1 
T ^ 22i6o 



(7) 



(8) 



From this result we see that the polymers become stiffer 
if one decreases the weights on (one type of) the non- 
vertical bonds, hence preventing transverse wandering. 



III. CLEAN SYSTEM: CONTINUUM LIMIT 
AND THERMODYNAMIC ENSEMBLES 

Before treating the random system, let us first com- 
pare the free energy density of the dimer models and the 
corresponding polymer lattice by taking the continuum 



3 



^ I ^ 




































1 














1 
































1 


















1 











































1 


1 




























1 




























1 


























1 
















FIG. 1: Mapping of dimer configurations to directed polymers 
by superposition with a fixed reference covering (middle). 



limit. The free energy of dimer model on the hexagonal 
lattice can be computed exactly By changing variables 
from zi, Z2 to ?7 = Zi/z2 and poj one can express the free 
energy density in terms of the physical quantities of the 
polymers on the latticeii. The result is 
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where the energy is measured relative to the line zi + 
Z2 = 1 of vanishing polymer density. For this expression 
we can take explicitly the continuum limit by sending 
bo ^ while keeping the polymer density po fixed, i.e., 
we adjust zi and Z2 so that the arcsin in Eq. jSjl tends to 
zero Bq while 77 = z\j z^ is kept fixed. This yields the 
free energy density for the continuum version of the dimer 
model, expressed in terms of the polymer parameters. 
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where we used Eq. JTJ. 

An alternative approach to treat a system of interact- 
ing polymers in (1 + 1) dimensions is to regard each poly- 
mer as a world line of a boson in imaginary time, leading 
to the quantum theory of interacting bosons in one spa- 
tial dimension. Choosing the imaginary time direction 
along the occupied bonds in the dimer reference state, 
the non-crossing constraint is naturally implemented by 
a hard core repulsion between the bosons which in one 
dimension is equivalent to non-interacting fermions. The 
Pauli principle then automatically prevents crossing. If 
the length Ly of the polymers tends to infinity, their re- 
duced free energy LyE^jh is given by ground state en- 
ergy Eq = [tt"^ /Q){'h? lm)p^Lx of ID fermions at density 
p. Using the mapping m ^ g and h ^ T, one gets for 
the reduced free energy density of the polymers at fixed 
density 
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Although the scaling of this result with the physical pa- 
rameters is the same as for the free energy in the contin- 
uum limit of the dimer model, Eq. (|l()(l , the amplitudes 
do not agree. However, since in the dimer model the 
number of polymers is not fixed, we have to compare the 
dimer free energy with the potential of the grand canon- 
ical ensemble of the polymers. The chemical potential 
is obtained as /i = dfpoiy/dp = {tt^ /2){T/g)p^ , yielding 
the grand canonical potential density 
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which is in full agreement with the continuum limit of 
the exact solution of the dimer model of Eq. 1)10(1 . This 
demonstrates that dimer model can be described on large 
length scales as free fermions where their mass is deter- 
mined in terms of the bond weights by a random walk of 
a single polymer on the lattice. 

It is instructive to compare the exact lattice result for 
the dimer free energy of Eq. ^ and the potential jpoiy 
with go of Eq. ((T)) even for larger bopo- By numerical in- 
tegration of Eq. we obtain the ratio /dimcr/jpoiy over 
the entire range of possible polymer densities shown in 
Fig. 121 Up to approximately 1 /4 of the maximal density, 
we find reasonable agreement between the lattice and 
continuum results, almost independent of the anisotropy 
r] — z\l Z2. For larger densities the value of 77 becomes im- 
portant. There is an optimal value of 77 close to 1/3 for 
which the continuum description gives accurate results 
(within a few percent) even for all densities. 

The anisotropy of the dimer model can be tuned by 
changing the relative magnitude of the weights z\ and z^- 
The exact solution of the dimer model on the hexagonal 
lattice yields for the correlation lengths the result-^® 
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with (f)o — arcsin {{zi + Z2Y' ^ l)/(4ziZ2)- Hence the 
correlation length perpendicular to the direction of the 
polymers is set by their mean distance 1/po- The length 
^j, should then be set by the typical scale a polymer can 
wander freely before it reaches a transverse displacement 
of the order of the mean distance between polymers. In 
the continuum description of the polymers, the random 
walk description of a single polymer then implies the re- 
lation 
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between the correlation lengths. Together with the first 
relation of Eq. (|13|l , this yields the anisotropy 
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That this result is consistent with the anisotropy of the 
elastic description of the polymer system follows from 
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FIG. 2; Comparison of the free energy /dimcr of the hexag- 
onal lattice dimer model, Eq. @, and the grand canonical 
potential jpoiy of the continuum polymer system, Eq. H12|l 
with the expression of Eq. 10 for g/T substituted. The 
curves extend to po = 2/(3^''^6o) which is the maximal 
polymer density on the lattice for rj = 1. However, in 
the lattice model, the density cannot become larger than 
PC = 2/(7r\/3) arcsin(-y2 + 77/2) by tuning the weights at a 
fixed ratio 77. 



model, we can use the insight we gained from the previ- 
ous analysis of the hexagonal lattice to compare the free 
energies of the square lattice and the continuum model. 
The exact reduced free energy density of the lattice model 
is known to bei 
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The continuum free fermion result of Eq. (|12() yields in 
combination with the random walk result of Eq. ((HJ for 
the line tension on the square lattice the estimate for the 
reduced free energy density 
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Below we will study the square lattice with bond ener- 
gies that are randomly distributed with mean zero so 
that the clean limit corresponds to the isotropic case 
zi = 1. For this case Eq. (|17|l yields the exact re- 
sult /dimor^n = — G/tt = —0.2916 which is close to 
the continuum approximation of Eq. p8|l which pre- 
dicts /dimcr^o = — 7r^/36 = —0.2742 at the fixed density 
Pn = l/(25n). 



IV. RANDOM BONDS AND PINNED 
POLYMERS 



the macroscopic compression modulus which is given by 
the compressibility, i.e., cn = T fpoiy / d in terms 
of the reduced free energy of Eq. (|ll|l . This yields cn = 
-k"^ {T"^ / g) p^ and hence together with C44 — gp 




In order to connect the lattice and continuum descrip- 
tion further, one would like to know under what condi- 
tions the lattice correlations lengths of Eq. H13|) fulfill the 
continuum relation of Eq. (|14|l . To address this question 
we change again variables from zi, Z2 to 77 = zi/ Z2 and 
PQ. If one uses the result for go of Eq. ([TJl one can 
easily check that Eq. H14() is indeed fulfilled in the con- 
tinuum limit 60 ^ 0. Hence, the exact anisotropy factor 
\/ C11/C44 of the continuum elastic model is recovered. 
This is important if one compares free energy densities 
of systems with different anisotropics since then the ratio 
of system sizes must be chosen as to match the ratio of 
their anisotropics. 

For the dimer model on the square lattice, there is 
no direct analog of the previous analysis since the mean 
polymer density cannot be tuned by changing the weight 
zi but is fixed, see Eq. ©. Hence, one cannot take the 
continuum limit explicitly. Nevertheless, the square lat- 
tice is particular useful if one wants to study the effect of 
random bonds since the polymer density is robust against 
variations of the weights and thus shows no disorder in- 
duced fluctuations. For the clean square lattice dimer 



In this section we will consider exclusively the square 
lattice but with random energies assigned to all verti- 
cal bonds so that Zij = exp(— ey/T^), where denotes 
now the weight on the bond {ij). The dimer tempera- 
ture measures the strength of disorder. The energies 
Cij are drawn for each bond independently from a Gaus- 
sian distribution with zero mean and unit variance. On 
all horizontal bonds we set e^j = so that for Td 00 the 
isotropic clean dimer model is recovered. The partition 
function can be written as 

= ^exp - J2 'v/Td] , (19) 
{D} \ iij)eD J 

where the second sum runs over all occupied bonds. 
The disorder averaged free energy and correlations 
of this model have been computed by a polynomial 
algorithmiLiS, for system sizes up to 512 x 512 lattice 
sites and typically 6000 disorder samples^iii. 

On the analytical side, progress has been made for 
the polymer system with random pinning by applying 
the replica method. Regarding again each polymer of 
the replicated theory as a fermion in imaginary time, 
and applying the Pauli principle for all particles within 
the same replica, the replica free energy can be obtained 
again as the ground state energy of a one-dimensional 
system of fermions. Due to the replication, the fermions 
carry now n spin components and interact via an attrac- 
tive (5-function potential arising from the short ranged 
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disorder correlations. This SU{n) fermi gas can be stud- 
ied by a series of nested Bethe Ansdtz^. In the hmit 
n — > 0, the Bethe Ansatz equations can be solved exactly 
for arbitrary disorder strength A, yielding the disordered 
averaged reduced free energy density of the polymersiS, 



/poly = /o(A)p 
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where /o(A) represents the disorder dependent free en- 
ergy of a single polymer. Notice the simple form of the 
disorder contribution to the free energy of the pure sys- 
tem, cf. Eq. Interestingly, in the limit of strong 
interactions (disorder) the SU{n) fermi gas in the limit 
n — !■ becomes identical to the (pure) interacting Bose 
gas studied by Lieb and Liniger^''. Since it was shown 
that the interaction strength scales as n^, perturbation 
theory for the ground state energy of the Bose gas yields 
a series expansion in n of the replica free energy for large 
disorder. The coefficients of this expansion correspond to 
the disorder averaged cumulants of the free energy which 
hence are known exactly from the replica Bethe Ansat2^. 

The prediction of the replica approach can be com- 
pared to the numerical evaluation of the free energy and 
its cumulant averages for the random bond dimer model. 
This has been done in Ref. |^ neglecting however fluctu- 
ations in the polymer density induced by the statistics of 
the pure dimer model. While nice agreement was found 
for the second and third cumulant of the free energy, the 
averaged free energies were only consistent if a term ~ A^ 
of the single polymer contribution /o was dropped by 
hand. However, contributions ~ n'^A^ from the replica 
free energy of a single polymer were found to be crucial 
for the agreement of the third cumulant of the total free 
energy. A similar observation was made^ for the data ob- 
tained previously for a single pinned polymer—. This is 
in particular unsatisfying due to the model character of 
the directed polymer in a random potential for the theory 
of disordered systems. Below, we show that the differ- 
ences between the replica approach for polymers and the 
numerical results on the dimer model can be fully rec- 
onciled when polymer density fluctuations are included. 
As demonstrated for the clean system, this can be done 
by comparing the canonical dimer ensemble to the grand 
canonical ensemble of polymers. From Eq. (|20|l follows 
the disorder averaged chemical potential pL = dfpoiy/dp 
since disorder induces no (additional) fluctuations in p. 
Thus, the average grand canonical potential density is 
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so that the single polymer term /o cancels. Notice that 
/o is exactly the term which was in disagreement with nu- 
merical results for the random bond dimer model. Hence, 
the disorder averaged free energy of a single directed 
polymer cannot be determined from numerical compu- 
tations of the free energy of the dimer model. How- 
ever, it is the disorder induced effective interaction of 



the polymers which determines the dimer free energy. 
To obtain the latter, we substitute the dimer parameters 
p ^ = l/(26n) and g/T from Eq. © into Eq. (EH- 
The disorder strength A/T^ must be related to the dimer 
temperature Td which measures the variance of the bond 
energies. As was shown in Ref. the relation is 
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where the length acts as a cutoff in the continuum 
model over which the 6 function of Eq. Q is smeared 
out. The ratio ^d/bn must be considered as a fitting pa- 
rameter which should turn out to be of order unity. As 
we are comparing energy densities, the disorder depen- 
dent anisotropy C11/C44 of the polymer system must 
be included^. This yields for the disorder averaged free 
energy density of the dimers 



/dimer — \l Jpoly • 



(23) 



The compression modulus can be obtained again from 
the (averaged) polymer free energy of Eq. ipn|) . cn — 
Tp'^d'^fpoiy/dp^, yielding for the anisotropy 
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In order to correct for a (small) difference between lattice 
and continuum model, we match with the exact result 
/dimer = — G/tt of Eq. (|17() with zi = 1 in the clean 
limit Tc( — > 00. Then we obtain in terms of the dimer 
parameters the final result 
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which has to be compared to the result for \nZi/L^ = 
— /dimer of Eq. (32) in Ref. 5. The expression of Eq. 
is plotted in Fig. |31 together with simulation data for the 
random dimer model, demonstrating indeed nice agree- 
ment for £,d/ba = 1.33. 

Finally, we comment on higher cumulant averages of 
the dimer free energy. They were also measured in sim- 
ulations, and were shown to agree with the free energy 
fluctuations of the polymer system^. This can be easily 
understood from the fact that on the square lattice there 
are no sample-to-sample variations of the mean polymer 
density. Hence the shift of the polymer free energy by 
—fip in Eq. H12|) is independent of disorder so that one has 
identical disorder averaged cumulants in the canonical 
and grand canonical ensembles, [j^]c = [f^]c for p > 2, 
where [. . .]c denotes a cumulant average over disorder. 
Because of that, fluctuations of the single polymer free 
energy /g are important for the dimer model. This ex- 
plains why contributions ~ n'^A^ from a single polymer 
to the replica free energy had to be included in Ref. 
to obtain agreement for the third cumulant of the free 
energy between polymer and dimer model. 
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FIG. 3; Disorder averaged free energy density for the ran- 
dom bond dimer model, Eq. 12511 with £,d/bn = 1.33, and 
corresponding simulation data, taken from Ref. 0. 



V. OUTLOOK 

We have shown that a continuum polymer model can 
provide a good approximation to dimer models on bipar- 
tite lattices. Although the polymer description yields in 
general not the exact result, it provides a more phys- 
ical picture of the dimer model as its exact solution 
which, moreover, is available only in the clean limit. Lat- 
tice Ising spin models with geometric frustration can be 
mapped at zero temperature exactly to classical dimer 
models on the dual latticeiS. The effect of thermal fluc- 
tuations and/or a transverse magnetic field can be un- 



derstood in terms of topological defects in the polymer 
representation of the dimer model^^ . Here we have shown 
that the influence of random bonds in the dimer model 
can be described as pinning of polymers. This implies 
that the glassy state of certain spin models with random 
couplings could be related to the glass phase of polymers 
in a random environment. For example, it can be eas- 
ily checked that random dilution of the triangular Ising 
antiferromagnet leads to a pinning of the polymers at 
the non-magnetic lattice sites. Another potential appli- 
cation of our results is the study of classical dimers which 
in addition to the hard-core repulsion interact in a more 
general way. Since the mapping to polymers is indepen- 
dent of the dimer interaction, one can use the analogy 
between directed polymers and world lines of bosons in 
imaginary time to explore dimer interactions in terms of 
interacting bosons in one dimension. Recently, the clas- 
sical limit (without kinetic term) of the quantum dimer 
model on the square lattice has been shown to have a 
phase transition between a critical and a columnar phase 
due to the aligning interactioniSi. The correlations in the 
critical phase are found to decay with an exponent that 
varies continuously with the interaction amplitude. Us- 
ing the mapping to world lines of bosons, that exponent is 
determined by the compressibility of the bose gas which 
presumably can be modeled by a tight-binding Hamil- 
tonian with an infinite on-site repulsion and a nearest 
neighbour interaction^^. 
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